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m . 

f"^ , We summarize previous works on the exact ground state and the elementary excitations of the 

exactly solvable BCS model in the canonical ensemble. The BCS model is solved by Richardson 
£SJ ' equations, and, in the large coupling limit, by Gaudin equations. The relationship between this two 

■i , kinds of solutions are used to classiffy the excitations. 

& : 

^ , I. INTRODUCTION 

<n : 

One of the most fundamental problems in a many body system is to know which are the ground state (GS) and the 
' (— i " low energy excited states, which determine the thermodynamic properties and the response to external fields. In most 
cases the excited states can be understood in terms of a collection of elementary excitations characterized by their 
statistics, discrete quantum numbers as spin, charge and dispersion relation. In the pairing model of superconductivity, 
proposed by Bardeen, Cooper and Schrieffer in 1957, this problem was solved long ago in the grand canonical ensemble 

O ( for a large number of particles 0,0- However recent studies, motivated by the fabrication of ultrasmall metallic 
grains, show that the grand canonical BCS solution deviates strongly from the exact numerical or analytical solution 
for systems with a fixed and small number of particles (for a review see Q)- Most of the previous studies have focused 
on the ground state of the BCS system, and some of its excitations. In this contribution we shall review further 

g . progress concerning the understanding of the full excitation spectrum Q . 

£ . II. THE GRAND CANONICAL BCS ANSATZ 

o : 

The BCS model of superconductivity is characterized by the energy levels of the electrons and the scattering 
potential between them. When the latter is a constant, the BCS model is exactly solvable a la Bethe [j| and 
integrable 6]. This is the model we shall used to study the GS and excitations. The Hamiltonian of the reduced BCS 
| is defined by Q 

Hbcs = \ ej a c\ a c ja - G^c} + cj 'y 'V • > (!) 

3,<r=± 3,3' 

O , 

where Cj t ± (resp. c| ± ) is an electron annihilation (resp. creation) operator in the time-reversed states \j, ±) with 
energies Ej/2, and G is the BCS dimensionful coupling constant. The sums in (JTJ run over a set of N doubly 
degenerate energy levels Ej/2 (j = 1, . . . , N). We adopt Gaudin's notation, according to which Ej denotes the energy 
of a pair occupying the level j 0. To fix ideas we shall use the so called equally space or picket fence model, 
that is the one employed in the study of ultrasmall superconducting grains Q. This model is given by the choice 
Ej = d(2j — N — 1), where d = ui/N is the single particle energy level spacing and ui/2 is the Debye energy. The 
W . coupling G can be written as G — gd, where g is dimensionless. 

The BCS ansatz for the GS of this Hamiltonian in the grand canonical ensemble is given by Q, Q 

> 

X 

where |0) is the Fock vacuum of the electron operators and Uj, Vj are the BCS variational parameters given by 



i 



\BCS)=H(uj + Vjcl + cl_)\0), (2) 



£. = e . -e -G, Ej = J$ + A2. (4) 



* Talk presented at the 6th International workshop on Conformal Field Theory and Integrable Models, Chernologka, Russia, Sept. 2002. 
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In these eqs. £o and A are twice the chemical potential and the BCS gap, which are found by solving the following 
equations: 

where M is the number of electron pairs. The most studied case in the literature corresponds to the half- filled 
situation, where the number of electrons, N e = 2M, equals the number of levels N In this case the solution of 
eqs. © in the large N limit is given by A = 10/ sinh(l/g) and Eq = 0. 

The excited states in the grand canonical ensemble can be obtained acting on the GS ansatz (J2J with the Bogoliubov 
operators 7j, cr (a = ±) 

7ji,<xi ■ ■ ■ lj n ,<?„ \BCS), 7j,± = Uj c j± T Vj c] T , (6) 

and have an energy \{Ej 1 + . . . + E, J2 ). Recall that in our conventions Ej is twice the quasiparticle energy, thus the 
factor 1/2 in the energy of the states iJBJ. 



III. PAIR-HOLE REPRESENTATION OF BCS MODEL 



Every energy level j = 1, . . . , N has four possible states given by 

|0) : empty, 



cj + cj_ |0) : doubly occupied 



cjjO) : singly occupied, (7) 

t J 



An important property of the Hamiltonian Q is the "blocking" of levels which are singly occupied, which means that 
these levels decouple from the rest of the system. This is a consequence of the equation 

Hbcs c]M = I ct CT + ct CT H B csW), (8) 

where the state \tp) does not contain the operators ct . Thus the singly occupied levels only contribute with their 
kinetic energy and one can study the dynamics on those levels which are either empty or doubly occupied. 

Let us call Hn.m the Hilbert space of states of M pairs distributed among N different energy levels. To describe 
7~(-n,m let us define the hard-core boson operators 

l\, 'V <7-- &] = ct i+ ct_, Nj=b\bj, (9) 
which satisfy the commutation relations 

[bj,b],] = Sj, f (1-2AT,-). (10) 
The Hamiltonian restricted to the Hilbert space Hn,m can then be written as 

Fbcs = 5>^, - G £ b]b r . (11) 



A basis of states of Hn,ai is given by 



i/>=n & ii°>' ( i2 ) 



where / denotes a set of M different integers ranging from 1 to N. Thus the dimension of Hn,m is given by_ the 
combinatorial number . A convenient pictorial representation of the singly and occupied states is given by 



ov 1 

a 

|0), .~&t| ). (13) 



At G = the ground state of (jl L|> is the Fermi sea obtained filling all the levels below the Fermi energy ep. The set Iq 
is given in this case by Jo = {1,2,..., M} (see fig. 1). Similarly the lowest energy excited state at G — corresponds 
to the choice I x = {1, 2, . . . , M - 1, M + 1} (see fig. 2). 
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GROUND STATE FOR N=8 LEVELS AND M=4 PAIRS 
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FIG. 1: Pair-hole representation of the ground state at G = for N = 2M 
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LOWEST EXCITED STATE FOR N=S LEVELS AND M=4 PAIRS 
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I,= {1, 2, 3, 5 



FIG. 2: Pair-hole representation of the lowest energy excited state at G = for N = 2M = 8. 

IV. EXACT SOLUTION OF THE BCS MODEL IN THE CANONICAL ENSEMBLE 

In 1963 Richardson showed that the eigenstates of the Hamiltonian with M pairs have the (unnormalized) 
product form [j| 



M N 

\M) = J] 4|vac), St=E -—jr b] , (14) 

8 7 H/j ■ 



where the parameters E u [y = 1, . . . , M) are, in general, complex solutions of the M coupled algebraic equations 

1 " 1 ™ 2 

G^J-E;- E E~E^ » = 1,-,M, (15) 

which play the role of Bethe ansatz equations for this problem 0, [tj IToL ITll IT2I IT3L ITiL Il5j . The energy of these states 
is given by the sum of the auxiliary parameters E v , i.e. 

M 

E{M) = ^E V . (16) 

The ground state of Hbcs is given by the solution of eqs. (|15|l which gives the lowest value of E(M). The (normalized) 
states (|14[) can also be written as [T^J 

l M > = 7=y E ^(ii J --. ) jM)6t i ...6t M | vac)) (17) 

' 31,— ,3 M 
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where the sum excludes double occupancy of pair states and the wave function ip takes the form 

M 

V fc=l 3k 

The sum in 118|l runs over all the permutations, V, of 1, • ■ • , M. The constant C in (|17|l guarantees the normalization 
of the state HU (i.e. (M\M) = 1). 

The number of solutions of the Richardson's eqs. I|15|) is equal to the dimension of the Hubert space TLn.Mi namely 
Cm 0- F° r finite and small values of N, M the solutions {E ll (G)}ff =1 have to be found numerically. These solutions 
and the corresponding eigenstates can be classified according to the values taken by E^G) at G = 0, namely 

Um^(G) = £j -, j el, (19) 

Or — ►U 

where / is the set / that label the basis (|12l) of Hm.m- This means that the spectrum of the BCS Hamiltonian follows 
an adiabatic evolution as a function of G [171 flS| For very small values of G all the roots EJG) of eq. \ib\ are real 
and close to their G = value given in eq. I|19|) . This happens for all the states labelled by /. 



Ground State solution 



In particular for the GS, i.e. Iq, as we increase G the real roots Em and Em-i, nearest to the Fermi level, approach 
from above and below the energy £m-i, and become equal to it at some critical value G = G c \. For G > G c \ the two 
roots Em and Em-i become a complex conjugate pair. Increasing G further one encounters a similar phenomena for 
the roots Em -2 and Em -3, and so on, until all the roots become complex (see fig. [J. 

In the limit where N is large, while MJN and u> = dN remain finite the complex roots E^ form an open arc T with 
end points £0 ± i A (see fig. 4) 0,0,0]. This gives an interesting geometrical meaning to the chemical potential 
(eo/2) and the BCS gap (A/2) for the exactly solvable BCS model. There are also roots which stay real staying in 
the segment (— w, £1) where £1 is the intersection point of T with the real axis. 
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FIG. 4: Plot of the roots _E" M for the equally spaced model in the complex plane [3. The discrete symbols denote the numerical 
values for M — 100. The continuous lines are the analytical curves obtained in |18|. All the energies are in units of uj. 

V. EXCITATIONS OF THE CANONICAL BCS MODEL 

In the canonical ensemble, where the number of electrons is fixed, the excited states can be obtained from the GS 
in two ways: by breaking Cooper pairs or by pair-hole excitations 0. 

Breaking Cooper pairs 

In fig. 5 we show at G = the excited state obtained by breaking the electron pair nearest to the Fermi level. The 
spin up electron remains in the same energy level, while the spin down goes one level up, producing two blocked levels. 
As shown in fig. 5 the problem becomes identical to that of two decoupled levels and a system with two less levels 
active for pairing interactions. Hence for G > the energy of this excitation is the sum of the single particle energy 
of the blocked levels plus the GS energy of the system with these levels being removed. This situation is similar to 
what happens when the system has an odd number of electrons or when it is under the effect of an external magnetic 
field 0. In both cases there are singly occupied levels, which only contribute with their non interacting energy. This 
sort of excitations are the ones that have been mainly studied in the literature. 



ff 
ff 
ff 



ff 
ff 

■ft- 



FIG. 5: Excited state obtained by breaking a Cooper pair at G = 0. On the rhs the singly occupied levels are blocked and 
decouple from the rest of the system. 
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Pair-hole excitations 



The pair-hole excitations consist in promoting a pair below the Fermi level to a level above it. At G = the 
resulting state can be viewed as a pair-hole excitation of the Fermi sea (see fig. 6). This sort of excitations belong to 
the Hilbert space Hn,m and were first consider in reference 0. It is clear that a general excitation consists in the 
combination of broken Cooper pairs and pair hole-excitations. The rest of this review will focuse on the latter ones. 



e« 

£ 7 

e 6 



H- 



e 4 - 



FIG. 6: Two possible excited states obtained by moving one and two pairs above the Fermi level. 

In fig. 7 we plot the real part of {E tJi (g)}f l L 1 for three states, labelled Joi-Th-^ an d a system with N = 40 energy 
levels and M = 20 pairs as a function of the BCS dimensionless coupling g from to 1.5. The states are given by the 
choices 

Jo = {1,..., 18, 19, 20}, h = {1,..., 18, 19, 21}, I 2 = {1, 18,20,22}. (20) 

The state Jo is the GS of the system and the pattern that follows ReE fl (g) is the same as in fig. 3. Notice that as 
g — > 00 all the roots become complex and escape to infinity. The state I\ is the lowest energy state of the system, 
and as we see in fig. 7, in the limit g — > 00, the root ^20(5) stays real and finite while the remaining M — 1 = 19 
roots escape to 00. Finally, the state I2 has two roots £"20(3) and E\g(g) that stay real and finite in the g — > 00 limit. 
This is a general feature of the numerical solutions of the Richardson's equations l|15fl . namely in the limit g — > 00 
there are Nq roots E fl (g) that remain finite, not necessarily real, while the M ~ Nq remaining ones go to 00. As we 
shall see below the number Nq of finite roots can be interpreted as the number of elementary excitations of a given 
excited state. 



Ground State (a) 1 excitation (b) 2 excitations (c) 

' 1 1 1 1 1 1 1 1 1 1 1 — _J 1 1 <~ 




FIG. 7: Real part of for the equally spaced model with M = N/2 = 20 pairs and Ng = 0, 1, 2 excitations Q- 
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FIG. 8: Excitation energies E exc = E - E G s < 14d for M = 20 pairs at half filling @|. There are 44 =13+26+5 states 
corresponding to Ng = 1,2,3 respectively. The particle-hole symmetry reduces these numbers to 25 = 7+15+3. 



Excitation energies 

To support this conjecture we have computed the excitation energy E exc = E — Eqs of some low lying excited 
states. Fig. 8 shows the energies E cxc for a system with M = N/2 = 20 as a function of g. At g — E cxc is simply 
given by the energy needed to lift the pairs from the Fermi sea up to the unoccupied energy levels, but as g increases 
the excitation energies quickly converge to asymptotas whose slope is given essentially by Nq, 

lim E cxc = Ng A, A - gw. (21) 

g^oc 

Moreover using the methods of Refs.0,0one can show in the large N limit that the excitation energy of a Richardson 
state is given by 

N G 

E exc = J2 V(E a - £o) 2 + A 2 , (22) 

where £o is twice the chemical potential, and the energies {E a } a 2 1 are the ones that remain finite in the g — > oo 
limit. In the latter limit one has A ~ goj and eq. (12211 becomes eq. 1)2 l[l. The excitation energy given by eq. 1)22(1 fits 
quite well the excitation energies of our prototype example (N = 40, M = 20), as shown in fig. HO 

In order to compare our results with the BCS standard solution let us consider the excitation energy of a real 

Cooper pair, hi the Bogoliubov approach, which is given by \J £j + A 2 (notice that A = 2Abcs)- The 

standard Bogoliubov quasiparticle with an energy | ^/e 2 + A 2 would have to be compared with excitations involving 

broken Cooper pairs. Since E a in eq. (1221 lies between two energy levels, with £j+\ — £j — 2d ~ 1/N (e.g. in fig. 03 
^20(00) = with e 2 o < £-20 < £21), E a = Ej + 0(1/N). Therefore, our theory is consistent within 0(1/N) corrections, 
as it is well known from the existing relation between the canonical and grand canonical ensembles. 

VI. COUNTING THE NUMBER OF STATES WITH N G ELEMENTARY EXCITATIONS 

We said above that the dimension of the Hilbert space TCn,m of the states with M pairs distributed into N 
different energy levels is given by the combinatorial number C^, which is also equal to the number of solutions of 
the Richardson's equations (|15|l . The next question is: how many of these solutions contain Ng roots that remain 
finite in the g — > 00 limit? Let us call cIn,m,n g that number which must obviously satisfy = J^n =0 ^n,m,n g - 
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In reference Gaudin made the following conjecture 

f C N 

dN,M,N G = S q N ° 



C% G ^ Q<N G <M 
N G > M 



(23) 



In table 1 we give some examples of this formula. Since do = 1, there is only one state where all the energies go 
to infinity as g — > oo, which is nothing but the GS of the system. 



N° of solutions 


dm 


dM-i 




dt = N - 1 


do = 1 


E^ finite 
Ep, infinite 


M 



M - 1 
1 




1 

M — 1 




M 



Table 1.- Classification of roots in the g — > oo limit. Here <J/v G stands for c?jv.m,jv g - 

This conjecture can be motivated as follows. Suppose that in the limit g — > oo there is a set {E a (oo)}^= i of N G roots 
that remain finite while the remaining M — N G ones escape to infinity. Then from eqs. 115|l one derives that the roots 
E a = E a (oo) satisfy 



JV 



° = E 



i 



3=1 



E n 



E 



E a -E n 



a 



1,...,N G . 



(24) 



These equations are due to Gaudin and they appear in the diagonalization of a spin chain model known as Gaudin 
magnets |2(j. The number of solutions of eqs. (|24|l was shown by Gaudin to be given by l|23l) 7]. In reference we 
checked numerically this Gaudin's conjecture for small systems. Furthermore we were able to show, in the large N 
limit, that for finite g the roots E a — E a (g) satisfy an "effective" Gaudin equations 



N 



° = E 



i 



N G 

E 



-RieMsj-Ev) ^ a R( Ep ){Ep~E a 



(25) 



with R(E) — ^{E — e ) 2 + A 2 . As g — > oo one has A ~ gu; and eqs. become eqs. i(2*Ij). 



VII. A CLASSIFICATION PROBLEM 



Gaudin's conjecture (|23|l gives the number dN,M,N G °f states with a given value of N G , but says nothing of how 
to find them. Recall that the solutions of Richardson's eqs. I|15|l are obtained by starting from a given initial state 
at g = 0, labelled by the set /, and then increasing the value of g. Hence the problem is to find for each state / 
which is the number of roots N G , which remains finite in the g — > oo limit. This defines the function N G (I) which, 
in physical terms, gives the number of elementary excitations present in the state. The total number of states / with 
N G = N G (I) is nothing but d N:M ,N G - 

Finding N G (I) is a highly non trivial problem because it connects the two extreme cases 5 = and g = oo. From 
a numerical point of view it is also a challenging problem due to the existence of several sorts of singularities in the 
merging and splitting of roots with no obvious a priori pattern, except for the GS. Despite of these difficulties it is 
possible to give a simple algorithm yielding N G (I), which has an interesting combinatorial interpretation in terms of 
Young diagrams [4|. Before we give this algorithm we need some preliminary definitions. 



Path representation of states and Young diagrams 

The first idea is to associate to every state J, in the Hilbert space Hn,m, a path 7/ in the lattice Z 2 . This is done, in 
the pair- hole representation of /, by associating to every empty level, o, a horizontal link in Z 2 , and to every occupied 
level, •, a vertical link, starting from the lowest energy state. Placing the initial point of 7/ at (0, 0), then its final point 
is at [M, N — M). The number of up and right going paths from (0, 0) and (TV — M, M) is given by C^, which equals 
the dimension of H-n,m- This means that the path representation is faithful. In fig. 9 we show the paths associated 
to the ground state Io and an excited state / at g = 0. The path jj is given by (0, 0) — * (0, M) — > (JV — M, M), while 
the path associated to any excited state I lies always below it. This fact makes possible to associate a Young diagram 
Yj to every state J, whose boundary is given by the paths 7/ and jj with their common links being removed (see 



FIG. 9: Path representation of the ground state To and an excited state /. Combining these two paths one gets the Young 
diagram Yj. The number of elementary excitations of the state I is given by Ng = 3, according to the rule 1261 1. 



fig. 9). The GS is associated to the empty Young diagram, i.e. Yj = 0. Using these definitions, the funtion Nq{I) 
is given by the following algorithm 0: 

Ng(I) = number of boxes on the longest SW — NE diagonal of Yj. (26) 

In the example shown in fig. 9 we get Ng — 3. The algorithm Ij26(l was proposed in reference [4j to explain a large 
amount of numerical simulations. As explained in Q the physical mechanism underlying eq. I|26[l is the collective 
behavior of the holes and pairs that occupy the closest levels to the Fermi level. There are also combinatorial reasons 
to support it (see below) , however an analytical proof is needed. 



(a) (b) 




FIG. 10: a) The path and Young diagram of an excited state with Ng = 3. The dotted line denotes the position of the Fermi 
level, b) Real part of 75 M . For g large enough there is a real root (1) and a complex root (2) in agreement with the result 
N G = 3. 

The non-triviality of i|26[l can be seen in fig. 10, which shows the evolution of the real part of the roots as a 
function of g. Only for sufficiently large values of <?, and after a complicated pattern of fusion and splitting of roots 
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(2 real roots <-> 1 complex root) one observes that, indeed, Nq = 3. It seems very difficult to explain this result on 
the basis of Bogoliubov's quasiparticle picture which only works properly for large number of particles in the grand 
canonical ensemble. 

We have yet no proof of eq. (|2(j|) . however there is the following consistency check based on Gaudin's conjecture 
(|23[l : cIn,m,n g must be equal to the number of Young diagrams, Yj, associated to the paths 7/ going from (0,0) to 
(M,N — M), which have Nq boxes on their longest SW-NE diagonal. The proof of this result uses the methods of 
Ref. 0. This result can also be proved using the RSOS counting formulas. (We thanks A. Berkovich for pointing out 
this fact). As an example we show in fig. 11 how the paths corresponding to the case N = 6, M — Nq — 2 can be 
mapped into RSOS paths. 

We would like finally to mention the close relationship between the path representation and the associated Young 
diagrams of the BCS states with the crystal basis used in the Fock space representation of the affinc quantum group 



U q (Sl(2)) |22j, |23j, |24| . At the moment of writting this paper we do not know if there exists a quantum group structure 
underlying the exactly solvable BCS model. This may indeed be the case given that the integrability of the BCS 
model can be explained from an inhomogenous XXZ spin chain with boundary operators 0, 0, 0, 0] . 



N=6, M=N G = 2 



1) ALL PATHS (0,0) — > (4,2) 



2) ADDING TWO BOXES 



(4,2) 



(0,0) 



3) PATHS WITH N G =2 4) MAPPING TO RSOS PATHS 




FIG. 11: 1) Shows the Cf = 9 paths from (0,0) to (4,2), 2) The paths corresponding to No = 2 are those below the shaded 
boxes, 3) We select only those paths with No = 2 and 4) after a rotation the Nq = 2 paths become RSOS paths on a Bratelli 
diagram. 



VIII. CONCLUSIONS 

We have shown in reference Q that the pair-hole excited states of the exactly solvable BCS model in the canonical 
ensemble can be interpreted in terms of elementary excitations with peculiar counting properties related to the Gaudin 
model, which have no counterpart in the BCS-Bogoliubov's picture of quasiparticlcs. 

The combinatorics involved in the counting of elementary excitations is similar to that of RSOS models, which 
suggests that these excitations are in fact solitons which were called "gaudinos" in Q . 

Some open problems are 1) the relation between these gaudinos and the Bogoliubov quasiparticles, 2) analytic 
proof of the algorithm for computing the number of gaudinos Nq(I) for each BCS state, 3) working out the physical 
consequences of these excitations in the canonical ensemble for small systems, 4) generalization of these results to 
higher spin representations and other Lie groups and supergroups. 
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